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Abstract —In this paper, we consider simultaneous wireless 
information and power transfer (SWIPT) in orthogonal fre¬ 
quency division multiple access (OFDMA) systems with the 
coexistence of information receivers (IRs) and energy receivers 
(ERs). The IRs are served with best-effort secrecy data and the 
ERs harvest energy with minimum required harvested power. To 
enhance physical-layer security and yet satisfy energy harvesting 
requirements, we introduce a new frequency-domain artificial 
noise based approach. We study the optimal resource allocation 
for the weighted sum secrecy rate maximization via transmit 
power and subcarrier allocation. The considered problem is non- 
convex, while we propose an efficient algorithm for solving it 
based on Lagrange duality method. Simulation results illustrate 
the effectiveness of the proposed algorithm as compared against 
other heuristic schemes. 

I. Introduction 

Orthogonal frequency division multiple access (OFDMA) 
gains its popularity due to its flexibility in resource allocation 
and robustness against multipath fading, and has become 
a promising multi-access technique for multiuser wireless 
communications networks. 

Recently, simultaneous wireless information and power 
transfer (SWIPT) has appeared as an appealing solution to 
prolong the lifetime of energy-constrained wireless nodes, 
by enabling them to receive energy and information from 
the same signal. SWIPT has drawn a great deal of research 
interests Q-El. For instance, in m, the authors studied two 
practical schemes for orthogonal frequency division multiplex¬ 
ing (OFDM) based SWIPT, namely, power splitting and time 
switching. With time switching applied at each receiver, the 
received signal is processed either for energy harvesting or 
for information decoding at any time. When power splitting is 
employed at the receiver, the signal is split into two streams, 
then processed for energy harvesting and information decod¬ 
ing, respectively. In El, the authors considered a multiuser 
OFDM system with some users to decode information, and 
the remaining users to harvest energy. 

Besides, due to the increasing importance of information 
security, substantial works have been dedicated to information- 
theoretic physical layer (PHY) security (e.g. a-Ei) , as a 
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Fig. 1. System model. 


complementary solution to the traditional cryptography based 
encryption applied in the upper layers. The authors in Q 
considered PHY security in an OFDMA system, with the goal 
of maximizing the rate of best-effort users subject to the secure 
data rate requirements of confidential users. In El, OFDM 
based wiretap channel was considered and the achievable 
secrecy rate with Gaussian inputs was studied. Artificial noise 
(AN) is another important method for improving PHY security 
by degrading eavesdroppers’ decoding capability il, El- In 
lISl . in order to assist secrecy information transmission, AN 
is transmitted into the null space of the channels of secrecy 
users to interfere with the eavesdroppers. In ||9l, the authors 
proposed a time-domain AN scheme by exploiting temporal 
degrees of freedom from the cyclic prefix in OFDM modulated 
signals, even with a single antenna at the transmitter. 

When PHY security is considered in SWIPT, AN not only 
enables secrecy information transmission but also becomes a 
new source for energy harvesting. There are only a handful of 
works that have studied the secrecy wireless information and 
power transfer by properly designing the beamforming vectors 
at the multi-antenna transmitter Qo), ini. Secrecy communi¬ 
cation in SWIPT over fading channels was also studied in 
(Ha. However, AN aided secrecy information transmission in 
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Fig. 2. Block diagram of an OFDMA transmitter with AN generation and the corresponding receiver with AN removal. 


OFDMA-based SWIFT systems has not been addressed in the 
literature, which motivates this worJc. 

In this study, we consider the resource allocation problem 
in the AN aided secure OFDMA systems with SWIFT, con¬ 
sisting of two types of receivers, i.e., information receivers 
(IRs) and energy receivers (ERs). We first propose a new 
frequency-domain AN method in OFDMA-based SWIFT to 
facilitate both secrecy information transmission and energy 
transfer to IRs and ERs, respectively. Specifically, independent 
AN signals are added over orthogonal subcarriers (SCs) and 
only the intended IR knows the AN transmitted over each 
corresponding SC and is able to remove it before decoding 
the information. Our goal is to maximize the weighted sum 
secrecy rate of the IRs subject to minimum harvested power 
requirements of individual ERs. The formulated problem is a 
mixed integer programming problem and thus is non-convex. 
We propose an efficient algorithm based on Lagrange duality 
method, which solves the problem optimally when the number 
of SCs becomes large. 

II. System Model and Froblem Formulation 
A. System Model 

We consider a downlink OFDMA-based SWIFT system 
with the new consideration of FHY security as shown in 
Fig. [U The system consists of one base station (BS) with a 
single antenna, K single-antenna receivers and N SCs. The 
set of receivers is denoted by /C = {l,...,Ar}, among which 
Ki receivers are IRs given by the set /Ci and the other K 2 
receivers are ERs given by the set /C2, i.e., /Ci U /C2 = /C. 
The set of SCs is denoted as J\f. Furthermore, we assume 
that for each IR, all other receivers (IRs or ERs) are assumed 
to be potential eavesdroppers. The BS is assumed to know 
the channel state information (CSI) of all receivers. This is 
practically valid since the IRs and ERs need to help the BS 
obtain their individual CSI for receiving required information 
and energy, respectively. 

We propose a frequency-domain AN generation and removal 
method in OFDMA systems, similar to that in ifT^ over the 
time domain. The scheme is illustrated in Fig. |2] and described 
as follows. A large ensemble of random Gaussian sequences 
are pre-stored at the BS, and the indices of the sequences 
are regarded as the keys. After SC allocation to IRs, the 
BS first randomly picks N sequences from the ensemble 


(each corresponds to one SC) and transmits each of their 
indices (keys) to the IR to whom the corresponding SC is 
assigned. As the random sequence (or AN) is only known 
to the intended IR but unknown to all the other receivers, 
any potential eavesdropper cannot have access to the random 
sequence used at each SC. Moreover, in order to prevent the 
eavesdropper from decoding the random sequence by long¬ 
term observation of the signal, the BS randomly picks new 
random sequences and transmits the corresponding keys in a 
secret manner to the desired IRs from time to time, using e.g. 
the method proposed in ifTsll . 

The transmit signal comprises of the transmitted data sym¬ 
bol Sk,n from the BS to IR /c on SC n and the AN bearing 
signal Zk,n for IR k, k G /Ci and n € Af. It is assumed that 
Sk,n and Zk,n are independent circularly symmetric complex 
Gaussian (CSCG) random variables with zero mean and unit 
variance, denoted by Sk,n ^ CJ\f{0, 1) and Zk,n ^ CJ\f{0, 1), 
which are also independent over n. 

The transmitted signal to IR fc at SC n is given by 

(1 kXk,n)Pk,n'^k,n “b Y ^k,nPk,n^k,n^ (1) 

where pk^n > 0 is the total power at SC n and 0 < ak,n < 1 
is the split power ratio to generate artificial noise to be added 
at SC n, with SC n assumed to be allocated to IR k. Notice 
that if pfe „ > 0 and ak,n = 1 for any SC n, then this SC is 
used only for energy transfer, i.e., there is no information sent 
over the SC. 

Let hk,n denote the complex channel coefficient from the 
BS to receiver k at SC n, and /3fc „ denote the eavesdropper’s 
complex channel coefficient. The downlink received signal at 
IR fe on SC n and that at a potential eavesdropper who is 
wiretapping IR k over SC n are respectively given by 

^k,n — hk^n^k,n “b ’t^k^ (2) 

dSk,n — l^k,n^k,n “b C/c, (3) 

where the noise Vk and Ck are assumed to be independent 
and identically distributed (i.i.d.) as CAf{0,(T^). Here, we let 
= maxfc/gAc.fe'/fc indicating that the consid¬ 

ered eavesdropper of receiver k is the receiver of the largest 
channel gain among all the other receivers on SC n. We 
assume that the OFDM symbols are time slotted so that the 
length of each time slot is comparable to the channel coherence 
time, i.e., the channel impulse response can be treated as 





















time invariant during each time slot. As a result, the BS can 
accurately estimate „ of all receivers and Pk^n on all SCs. 

With the aforementioned scheme, the AN can be removed 
at the desired IR at each SC but not possibly at any of the 
potential eavesdroppers. From (l)-(3), the received signals at 
IR k after AN cancelation and the “best” eavesdropper on SC 
n are further expressed as 

(1 Oi/e,n)p/c,n'5fc,n “t” (4) 

Ek,n — Pk,n^J (1 kXk,n)Pk,n^k,n “t” Pk.ny/^k,nPk,n^k,n “f ^k- 

(5) 


Here we can write the achievable information rate of IR k on 
SC n, which is given by 

(1 I ^/c,n I PkjiL j 


rk,n = log 2 1 + 


The decodable rate of the “best” eavesdropper on SC n is 
given by 

(1 |/5/i;,n I Pk, 


= l0g2 1 + 


+ ak,n\Pk,n\'^Pk,r. 


(7) 


The achievable secrecy rate for IR k on SC n is thus given 
by a 

K.n - rlX 

(1 - ak,n)\hk,n\'^Pk,r. 


I0g2 1 + 


- I0g2 1 + 


(1 Oik^n)\Pk,n\ Pk,r, 

ak,n\Pk,n\'^Pk,n + CT^ 

for all /c G /Cl and n G Af, where [•]“*■ = max(0, •). 
Corollary 1: 7?^ „ in ® can be further expressed as 


( 8 ) 
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(9) 


( 10 ) 


Proof: Please refer to Appendix A. ■ 

The weighted sum (secrecy) rate of all Ki IRs is given by 


T^sutn ~ ^ ^ '^k ^ ^ Xk,nRk,r. 
k^Ki nGJ\f 


( 11 ) 


where Wk is the positive weight of IR k, and Xk,n is the binary 
SC allocation variable with ccfc „ = I representing SC n is 
allocated to IR k and Xk^n = 0 otherwise. 

Next, for the ERs, the harvested power at each ER fc G /C 2 
is given by Q 

Qk — Ck ^ ^ f ^ ^ ^k,nPk,n | |/^/c,n| i (17) 

rtGA^ \k&Ki / 


'Note that for the case of n = 0^ we set [Xf. > +oo \f\hf. nP > 

l/3fc,nP and [Xk,n]'^ = 0 if \hk,-n.? < \lik,nP- 


where 0 < < 1 denotes the energy harvesting efficiency. 

Note that in the considered system, the ERs can harvest energy 
from all SCs while the IRs need orthogonal SC assignment for 
avoiding multiuser interference. 


B. Problem Formulation 

Our goal is to maximize the weighted sum rate of the IRs 
by optimizing transmit power and SC allocation as well as 
power splitting ratio at each SC, subject to the harvested power 
constraints of all ERs. The problem can be mathematically 
formulated as 


max^i?|„^ 

(13) 

Qk ^ Qk^k € /C 2 : 

(14) 

^ ^ ^ ^ Pk,n^k,n ^ -^max 
k^Ki 

(15) 

0 ^ Pk,n ^ -fpeak; VtZ G Af , k G /Ci 

(16) 

0 < o:k,n < 1, Vn € A/*, k G JCl 

(17) 

Xk,n € {0, l},VneA/',A:e/Ci 

(18) 

X ^k,n < l,Vn G Af, 

kGK.1 

(19) 


where P = {pk,n\ denotes the power allocation over SCs, 
X = {xfc „} denotes the SC allocation policy, and a = 
{ttfe „} denotes the power splitting ratio over SCs. In (fT4l) . 
Qk denotes the harvested power constraint for ER fc G /C 2 . 
In (fTSl l and (fTbl) . Pmax and Ppeak represent the total power 
constraint over all SCs and the peak power constraint over 
each SC, respectively. Einally, and (SSll constrain that 
any SC can only be assigned to at most one IR. 


III. Proposed Solution 

The problem (fTSl l is a mixed integer programming problem 
and thus is NP-hard. As shown in m, the duality gap 
becomes zero in OEDM-based resource allocation problems 
that include the problem (fOl l as the number of SCs goes to 
infinity due to the time-sharing condition. This implies that 
problem (fOl l can be solved by the Lagrange duality method, 
as will be shown in this section. Eirst, the Lagrangian function 
of problem (O is given by 

C(P,a,X,\,y) 

= ^ Xk,nRl,n - 7 ( X] X! ^k,nPk,n 

fcG/Ci n^J\f \kG]Ci uGM 

+ X! ^k{Qk — Qk) 

fee/Ca 

= X X “ 7 X X ^k,nPk,n 

fcG/Ci n^J\f fcG/Ci 

+ X X ^k,nPk,n I ^ ^ ^kCk\^k,n\ 

nGAf KkGKi / k^K-2 

^ ^ ^kQk 
fcG/Ca 















where 


where [Ai, A 2 ,Ai^^] and 7 are the Lagrange multipliers 
(dual variables) corresponding to the minimum required har¬ 
vested power constraints and the total transmit power con¬ 
straint, respectively. 

We then define V for given X as the set of all possible 
power allocations of P that satisfy 0 < pk^n < ^peak for 
Xk,n = 1 and pk^n = 0 when Xk,n = 0, 5 as the set of all 
possible X that satisfy constraints (fTSl l and (fT^ . and A as the 
set of all feasible cx that satisfy (Ell. Then, we can obtain the 
dual function for the problem (fOl i as 

p(A, 7 )= max £ (P, a, X, A, 7). (21) 

Pev{x),aeA,xeS 

The dual problem is given by 

min p(A, 7 ). (22) 

A^0,7>0 

From (l20l i. we can observe that the maximization in (l2^ 
can be decomposed into N independent subproblems. Hence, 
we can rewrite the Lagrangian as 

£ (P, a, X, A, 7) = ^ £„ (P„, ^ XkQk 

nGAf k£K,i 

+ 7-Pmax, (23) 

where 

£n {Pni ^n) — '^kPk,n ^Pk,n “t" Pk,n ^ ^ XkCk\^k,n\ • 

fce/Ca 

(24) 

And the subproblem for SC n is given by 

max £^ (P^, (Ayj, XA) • (25) 

Pn&vix),a„eA,x„es 

A. Joint Optimization for Power Allocation and Power Split¬ 
ting Ratio 

It is difficult to directly express the partial derivative of R^. ^ 
given in (Ell with respect to either „ or ak,n- However, as 
we have discussed in Corollary 1, ^ = 0 when 0 < pk,n < 

and „ = rk,„ - when pk,n > [A’fc,„]+. In 

each case, „ is differentiable with respect to pk^n and ak,n- 
Hence, we first find the optimal power allocation ^ and 
optimal transmit power splitting ratio ^ in both cases. Then, 
we select {pi ’^^at achieves the largest Lagrangian. 

]) The case of pk^n > [Xk,n\^: By using Karush-Kuhn- 
Tucker (KKT) conditions and combining it to the constraint 
(E), we can obtain the optimal al „ with given p^ „ as 


ai = \TL2\hk,n\^{al^^ - ak,n)\lik,nt^n, (28) 

~{^k,n tr/j; yj) yj I Wk 

+ ln2|/3fc ,[{al „ - - \Pk,n\‘^ak,n] 

(29) 

Cl =ln2(afc,„ - 

-b - ak,n)\Pk,n\'^\hk,n\'^WkCr‘^, (30) 

di ={ak,n - l){\hk,n\^ - |/3fc,„nw/ccr^ - ln2cr®H„, (31) 

= — 7 + ^ XkC,k\hk,n\^ ■ (32) 

k(ZK2 

We define $ 1 ( 0 ;^ „) as the set of all non-negative real roots 
to (El that satisfy [Xk,n\^ < Pk,n < Ppeak with given ak,n- 
To jointly optimize p^ „ and ak,n, we consider the following 
two subcases. 

For the first subcase, we remove the [-Jg operator and 
assume that al ^{pk,n) lies in [0,1]. Substituting it into (l27l i 
to eliminate ak,n, we have 

+ hPk.n + C2 = 0, (33) 

where 

02 =ln2|/3fc,„|4|VnPf2„, (34) 

&2 =tUfc|/3fc,nhfifc,np+ln2H„|/3fe,(35) 
C2 =0-^ 

+ ^"^‘^^l{\f3k,n\'^ + \hk,n\'^)} ■ (36) 

Similarly, we define $2 as the set of all non-negative real roots 
to (l3^ (also the feasible candidates) that satisfy [Afe „]+ < 

Pk,n if Ppeak- 

For the second subcase that al ^{pk,n) is greater than 1 or 
smaller than 0, the set of the feasible candidates for pi ^ is 
given by <hi(afc_„ = 1) U ^i{ak,n = 0), obtained via (El l. 

2) The case ofO < pk^n < [Xk,n\^: As we have discussed, 
„ = 0 in this case. The Lagrangian function can be thus 
rewritten as 

{Pn^ ^n) — Pk,n ^ ^ Afc | | ^Pk^ni (37) 

keK2 

which is a linear function with respect to pk,n and independent 
of ak,n- The feasible candidate pk,n in this case can be thus 
obtained as 

0, if 7 > Efcex;2 

min{[T’fc_„]+, Ppeak}, otherwise. 

(38) 


Pk,n{Pk,n) 


'1 (|/ 3 fc,nP - |fefc,nP)g^ ' ^ 

2 2\l3k,n\'^\hk,n\'^Pk,n . g ’ 


for all fc G /Cl and n G AC, where = min{max{-, a}, b}. 

On the other hand, by deriving the partial derivative of £„ 
with respect to pk,n and equating it to zero, we have 


+ blPk,n + ClPk,n + di = 0, 


Finally, combining the above discussions, we define 

T = ^l{ak,n = 0) U ^l{ak,n = 1) U $2 U {0, [<Tfc_„] + , Ppeak} 

(39) 

as the set of all feasible candidates for optimal pi The 
optimal power allocation pi ^ can be obtained as 

Pl,n = arg max £„ {pk,n,atn{Pk,n)) , 

Pk,n^J^ 


(27) 


( 40 ) 





Algorithm 1 Proposed Solution for Problem ( fT3] ) 

1: repeat 

2 : Find the set of all feasible candidates T according to 

(O and (l39ll. 

3: Compute the corresponding a\^{pk^n) for all candi¬ 

dates. 

4: Select the optimal and al ^iPk n) according to 

gii. 

5: Solve sc allocation ^ for all k € ICi and n € J\f 

according to (EB. 

6 : Update A and 7 according to (|4^ and (l44li respectively. 

7: until A and 7 converges. 


where ^{pk,n) is the optimal power splitting ratio corre¬ 
sponding to each optimal power candidate. 



Energy Harvesting Requirement (mW) 


B. Subcarrier Allocation 


Fig. 3. Achievable secrecy sum-rate -Rfum versus harvested power require¬ 
ment. 


Substituting the optimal ^ and ^ into the optimal 
SC assignment policy is given by 


^k,n 

where 


1 , if fc = fc* = argmax^gKi '^kAPk,n^o*k,n) 
0 , otherwise, 


(41) 


'kLk,n{,Pk,ny ^k,n) — '^kRk,n ^Pk,n “t” Pk,n ^ ^ • 

fcGKa 

(42) 


C. Dual Update 

According to ca, the dual problem is always convex; 
hence, the subgradient method can be used to update the dual 
variables to the optimal ones by an iterative procedure: 


= [^I- - {Qk - Qk)]'^ yk e /C2, 


T — T ^ f -^max ^ ^ ^ ^ ^k,nPk,n 

\ nGAf kGJCi / 


(43) 

(44) 


where f > 0 is the iteration index, and [Ci; ^ 

properly designed positive step-sizes. 


D. Complexity 

The complexity of this iterative algorithm is analyzed as 
follows. For each SC, 0{Ki) computations are needed for 
searching the best IR. Since the optimization is independent at 
each SC, the complexity is 0{KiN) for each iteration. Last, 
the complexity of subgradient based updates is polynomial 
in the number of dual variables K 2 + 1 Ea. As a result, 
the overall complexity of the proposed algorithm for solving 
problem (O is 0 {{K 2 + l)'^KiN), where q is a positive 
constant. 

Finally, we summarize the algorithm in Algorithm 1. 


IV. Numerical Results 

In this section, we evaluate the performance of the proposed 
algorithm through simulations. In the simulation setup, a single 
cell with radius of 200 m is considered. The BS is located at 
the center of the cell. The number of the SCs is TV = 64. We 
assume the noise power of = —83 dBm. We consider Ki = 
4 IRs that are randomly located in the cell with distance to the 
BS uniformly distributed. For each IR, rcfc = 1, V/c € /Ci. We 
also consider K 2 = 4 ERs that are uniformly distributed with 
radius of 2 m around the BS0 We also assume that the ERs 
all have the same harvested power requirement, i.e., Qk = 

0 ,Vfce/C 2 . 

Eor performance comparison, we also consider several 
benchmarking schemes. Eirst, the fixed transmit power split¬ 
ting ratio with a = 0.5 or a = 0.2, is considered for 
complexity reduction, where the power and SC allocation is 
still optimized as the proposed algorithm. Here we drop the 
index k, n in ak,n for brevity. Second, the SC assignment 
is fixed (ESA) while the power allocation and transmit power 
splitting are jointly optimized as the proposed algorithm. Last, 
we consider the scheme without using AN (NoAN). 

In Eig. [3] the sum rate versus harvested power 

requirement Q is demonstrated with Pmax = 37 dBm. Eirst, 
for all schemes (except NoAN), the sum rate decreases with 
increasing Q (except NoAN). An interesting observation is 
that the scheme with a = 0.5 performs closely to the proposed 
algorithm. Moreover, the proposed schemes with AN achieve 
great rate-energy gains over that of NoAN, which has almost 
zero secrecy rate regardless of harvested power requirement. 
This is because without AN, the secrecy rate on each SC is 
positive only when it is assigned to the receiver of largest 

^We consider ERs in general closer to the BS than IRs to receiver larger 
power (versus that of IRs used for decoding information against background 
noise only). However, under this circumstance, the ERs have better channel 
conditions than the IRs, and as a result they are more capable of eavesdropping 
the information □S). 
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Fig. 4. Achievable secrecy sum-rate Rsum versus total transmit power Pmax- 


one root at „ = 0, when > \Pk,n\'^, and two roots 

at pk,n = 0 and pk,n = Xk,n, when \Pk,n? > 

Furthermore, when it follows that 


d{rk,n - rl ) 


dpk,n 


\Pk,Ti=‘^k,n 

ak,ni\Pk,n\'^X^ri + ^k,nCr^) 


> 0 . 


Xk.r^ 

(46) 


Hence, rk,n - < 0 when 0 < pk,n < Xk,n and 

\Pk,n? > \hk,n?, which is equivalent to 0 < pk,n < [Xk,n\^■ 
'^k,n - > 0 when Pk,n > Xk,n, if 1/3/c. np > |3l/c.nP 

or pk,n > 0 if > \Pk,n\'^, which is equivalent to 

Pk,n > [Xk,n\^■ We finally show that ^ = 0 when 
0 < Pk,n < while = rk,n - > 0 when 

Pk,n ^ [3t/c^7^]~*". 

The proof is thus completed. 


channel gain Q. In our simulation setup, the ERs possess 
much better channel gains compared to the IRs, due to much 
shorter distances to the BS. As a result, P < l/3/c.nP is 
true for all n € Af, fc € /Ci in general, and hence no secrecy 
information can be transmitted at all. This demonstrates the 
effectiveness of the proposed AN based approach. 

Fig. El demonstrates the sum rate R^nm versus the total 
transmit power Pmax. with the harvested power requirements 
sets as (5 = 100 /rW. Compared with the other benchmarking 
schemes, the proposed algorithm is observed to perform the 
best. 

V. Conclusion 

This paper studies the optimal resource allocation problem 
in OFDMA-based SWIFT with the new consideration of PHY 
security. With a proposed frequency-domain AN generation 
and removal method, we maximize the weighted sum rate 
for the secrecy IRs subject to individual harvested power 
constraints of ERs by jointly optimizing transmit power and 
SC allocation as well as transmit power splitting ratios over 
SCs. The formulated problem is solved efficiently by a pro¬ 
posed algorithm based on Fagrange duality. Through extensive 
simulations, we show that the proposed algorithm outperforms 
other heuristically designed schemes with or without using the 
AN. 

Appendix A 
Proof of Corollary 1 

Equating rk,n — ft n to zero, we have 

\hk,n\^Pk,n |/3/i;,n| Pk,n ('4^1 

a/c.n|/3/c,„pPfe,„ + Cr2' 

The roots of the above equation are given by „ = 0 and 
Pk,n = = 'Va:,™- However, since pk,n is 

always non-negative, pk,n = 44 „ > 0 can be true only when 
|/3/£.nP > \ hk , n \'^- Thus, we show that rk,n - = 0 has 
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